#### lecture 13, Fig. 5 ######
curve(dnorm,-4,4)
xvar<-seq(-4,2,length=100)
dvals<-dnorm(xvar)
## pnorm
polygon(c(xvar,rev(xvar)),c(rep(0,100),rev(dvals)),col='grey70')
arrows(1,.15,2,.25,code=1,angle=45,length=.1,lwd=2)
text(2,.25,paste('pnorm(2) = ',round(pnorm(2),3)),pos=3,cex=.85)
### dnorm
segments(2,0,2,dnorm(2),lwd=2,col='tomato')
arrows(2,dnorm(2),2.35,.13, code=1, angle=30, length=.1, lwd=2, col=2)
text(2.35,.13,paste('dnorm(2) =', round(dnorm(2),3)), cex=.85, pos=3, col=2)
### qnorm
points(2,0,col='dodgerblue',pch=16,cex=1.1)
arrows(2,0,3.1,.05,code=1,angle=30,length=.1,lwd=2,col=4)
text(3.1,.05,paste('qnorm(',round(pnorm(2),3),') = 2'), cex=.85, pos=3, col=4)
#rnorm
par(xpd=T)
points(rnorm(50),jitter(rep(0,50)),pch=16,cex=.6,col='seagreen')
arrows(-1,0.01,-2,.15,code=1,angle=30,length=.1,lwd=2,col='seagreen')
text(-2,.15,'rnorm(50)',col='darkgreen',cex=.85,pos=3)
par(xpd=F)
mtext(side=3,line=.3,'X ~ Normal(0,1)',cex=.9,font=2)